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The magnetic behavior of insulating doped diluted magnetic semiconductors (DMS) is character- 
ized by the interaction of large collective spins known as bound magnetic polarons. Experimental 
measurements of the susceptibility of these materials have suggested that the polaron-polaron in- 
teraction is ferromagnetic, in contrast to the antiferromagnetic carrier-carrier interactions that are 
characteristic of nonmagnetic semiconductors. To explain this behavior, a model has been devel- 
oped in which polarons interact via both the standard direct carrier-carrier exchange interaction 
(due to virtual carrier hopping) and an indirect carrier-ion-carrier exchange interaction (due to the 
interactions of polarons with magnetic ions in an interstitial region). Using a variational proce- 
dure, the optimal values of the model parameters were determined as a function of temperature. 
At temperatures of interest, the parameters describing polaron-polaron interactions were found to 
be nearly temperature-independent. For reasonable values of these constant parameters, we find 
that indirect ferromagnetic interactions can dominate the direct antiferromagnetic interactions and 
cause the polarons to align. This result supports the experimental evidence for ferromagnetism in 
insulating doped DMS. 



PACS numbers: 75.50.Pp, 75.10.-b, 75.30.Hx 
I. INTRODUCTION 

Diluted magnetic semiconductors (DMS) are semicon- 
ductors in which a fraction of the nonmagnetic ions that 
make up the crystal structure have been replaced by 
magnetic transition metal or rare earth ions. For ex- 
ample, substituting Mn^+ ions for some of the Cd ions 
in the nonmagnetic semiconductor, CdTe, yields the di- 
luted magnetic semiconductor, Cdi_a;Mn2,Te. In doped 
DMS, the sizable exchange interaction between magnetic 
ions and carriers (electrons or holes) leads to unusual op- 
tical, magnetic, and transport properties. Due to their 
potential for use in novel devices which take advantage 
of both their semiconducting and magnetic properties. 
DMS have, of late, beeit-the subject of much interestBEl 
Recently, the discoveryEl of a ferromagnetic transition 
temperature of 110 K in a sample of Gai_a;Mn2;As with. 
X « 0.05 has fustier enhanced both the experimentaloEl 
and theoretical Il3 interest in DMS. 

In the II-VI DMS Zni.^^Mn^Te {x < 0.1), p-doped 
with carriers at the level of 3 x 10"'^ "^/cm"^, where the sys- 
tem is in the insulating state, measurements of suscep- 
tibility versus applied magnetic field were conducted by 
J. Liu at NEC Research Institute. Tie data, originally 
reported in an NEC technical memolij, has been repro- 
duced in Fig. 1^ for the convenience of the reader. Note 
the double-step structure of the susceptibility and the 
two characteristic field scales indicated by the inflection 
points of the curve. This form suggests a dual magne- 
tization mechanism whereby large collective spins align 
at fields 300 G) too weak to magnetize the individual 
magnetic ions. Only at much larger fields {^-^ 15,000 G) 
do the individual Mn spins align. Within this interpre- 
tation, the measured susceptibility can be viewed as the 



sum of two contributions: a collective spin term that 
drops off around 300 G and an individual spin term that 
drops off around 15,000 G. The dashed Hne in Fig. | 
serves to separate these two contributions. 

For the x-regime in question, the Mn concentration is 
not enough to percolate, and the undoped system is not 
magnetically ordered. (Spin glass type order of the un- 
doped system has been obser-jiad in II-VI DMS for Mn 
concentrations above x — 0.2.113) Consequently, the un- 
usual magnetic behavior is attributable to the presence of 
the dopants. This is not surprising, despite the low car- 
rier concentration, because the Bohr radius which char- 
acterizes the carrier wave function is large compared to 
the Mn 3d wave function which characterizes the extent 
of the Mn local moment. We interpret the large collec- 
tive spins, responsible for the double-step form of the 
susceptibility, to be bound magnetic polarons, formed by 
the exchange interaction between localized carriers and 
magnetic ions within the carrier orbit. Furthermore, a 
fit of the polaron part of the the susceptibility data to a 
Curie- Weiss form reveals a net ferromagnetic interaction 
between the polarons. This result is in stark contrast 
to that observed for conventional nonmagnetic semicon- 
ductors in which virjtaal carrier hopping invariably yields 
antiferromagnetismJla To explain both the formation of 
bound magnetic polarons and the ferromagnetic nature 
of their interaction, we introduced, in Ref. a bound 
magnetic polaron model for insulating doped DMS. This 
was further elucidated by a comprehensive calculation in 
which we showed how the parameters of the model could 
be obtaipfid in an optimal manner using a variational 
principles, which we present below. 

In Sec. |l| we describe the system of two interacting po- 
larons in a diluted magnetic semiconductor and develop 
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FIG. 1. Magnetic susceptibility [dM/dH) of Zni_a;Mn^Te 
{x « 0.1), p-doped at the level of 3 x 10^^/cm^, measured at 
2 K as a function of applied magnetic field. Circles denote 
measured data while the solid line is a fit to a dual mag- 
netization model. The dashed line separates the collective 
spin (polaron) contribution from the indiwdual Mn spin con- 
tribution. Data was obtained by J. LirO at NEC Research 
Institute and is reproduced here with the permission of NEC. 



a simplified model to capture its behavior. In Sec. 
we calculate, within our model, both the single polaron 
partition function and the polaron-pair partition func- 
tion. Making use of these partition functions, we im- 
plement a variational procedure, in Sec. IV, to optimize 
the parameters of our model. We find that while the 
model parameters describing single polaron formation are 
temperature-dependent, the polaron-polaron interaction 
parameters can be treated as temperature-independent 
constants for magnetic ion and carrier densities of in- 
terest. We make use of these results in Sec. ^ where 
we demonstrate how a ferromagnetic polaron-polaron in- 
teraction can come about. Conclusions are presented in 
Sec. IVI. 



II. POLARON-PAIR SYSTEM AND MODEL 

A. The System 

To understand the magnetic behavior of diluted mag- 
netic semiconductors, we consider the polaron-pair sys- 
tem which consists of two carriers (electrons or holes) 
bound to impurity sites (donors or acceptors) separated 
by an inter- impurity distance, -R12, and the magnetic 
ions (usually spin 5/2 Mn) which surround them. This 
complex system interacts via three independent exchange 
interactions each of which have different characteristic 
length scales. 



The bound carriers interact directly via an impurity- 
impurity exchange interaction. Although this interaction 
can be more complicated for the casCpOf acceptors (for 
which the valence band is degenerate)Lj, we assume an 
impurity-impurity interaction of the Heisenberg type, as 
in the donor case, characterized by an exchange constant 
J. This interaction has been shown to be antifeiromag- 
netic for donors in nonmagnetic semiconductorsEEl and is 
assumed to be so for carriers in diluted magnetic semi- 
conductors as well. J is a function of the inter-impurity 
distance and the effective Bohr radius and varies as 



J exp(-2i?i2/as). 



(2.1) 



Were it not for the influence of additional interactions, 
this direct exchange would yield a net antiferromagnetic 
exchange interaction in DMS. 

The second interaction at work in this system is the 
exchange interaction between each of the carriers and the 
magnetic ions. This interaction is also antiferromagnetic 
and is proportional in magnitude to aj^'p where a is the 
carrier-ion exchange constant for the particular material 
and ^E* is the carrier wave function. For the purpose of 
this study, we will take 4' to be the hydrogenic wave 
function 



^'(r) = (7ra|)-i/2exp(-r/as) 



(2.2) 



with an effective Bohr radius a^. However, it should be 
noted that for acceptors in particular, the carrier wave 
functions may be more complicated.Ej We are also im- 
plicitly assuming that the binding energy of the impurity 
is large compared to the magnetic energy of the polaron, 
so magnetic ordering does not change the carrier wave 
function. 

Finally, there exists an additional antiferromagnetic 
exchange interaction, between the individual magnetic 
ions, which has a characteristic length scale on the or- 
der of a magnetic ion radius (~ A). Since this length 
scale is small compared with others in the system (i.e. 
ttB ~ 10 — 20A), we neglect all but nearest neighbor in- 
teractions and assume that the nearest neighbors form 
inert singlets. Thus, ion- ion interactions are considered 
only via the use of an effective magnetic ion concentra- 
tion, X = x{l Y^)^^ , in place of x, the true magnetic ion 
concentration.im 

Hence, the polaron-pair system (given the assump- 
tions noted above) interacts via two antiferromagnetic 
exchange interactions, a carrier-carrier interaction and a 
carrier-magnetic ion interaction. These interactions are 
depicted in Fig. ^ and result in the Hamiltonian 

= a^si • S^l^-inp + a^S2 • S„|4'2„P + Jsi • S2 

n n 

(2.3) 

where n runs over all magnetic ions, and are the 
carrier wave functions at the magnetic site n, si and S2 
are the carrier spins, and the Sn are magnetic ion spins. 
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FIG. 3. Schematic of polaron-pair model. 



FIG. 2. Schematic of polaron-pair system. 



We consider this to be a polaron-pair system because 
carrier-ion interactions tend to anti-align the spins of 
magnetic ions in the vicinity of a carrier with respect 
to the carrier spin. Thus, each carrier and the ions in 
its vicinity form a single magnetic polaron with a large 
collective spin. The polarons interact via both the di- 
rect antiferromagnetic carrier-carrier exchange interac- 
tion and the indirect ferromagnetic exchange interaction 
that results when carrier-ion interactions cause both po- 
larons to anti-align with the same magnetic ions. Details 
of the competition between these two interactions will be 
explored as we study the nature of a simplified model. 



B. The Model 

Although the true Hamiltonian provides the best de- 
scription of the polaron-pair system, its solution is com- 
plicated by the fact that the magnitude of the carrier-ion 
interaction varies exponentially with carrier-ion distance 
(since carrier wave functions are hydrogenic). In order 
to obtain a more detailed understanding of the polaron- 
pair system, it is necessary to study a simplified, more 
tractable model. Such a model is obtained by making 
two simplifying approximations: the single step approxi- 
mation and the interstitial region approximation. 

The single step approximation entails replacing the 
carrier wave functions by radial step functions that are 
constant up to a radius R and zero beyond R. In this ap- 
proximation, all of the magnetic ions within a sphere of 
radius R about a carrier interact with that carrier with 
the same exchange constant, K. Thus, in this model, 
the definition of a polaron becomes clear. A polaron is 
composed of a single carrier and all of the magnetic ions 
within a radius R of the impurity site to which the car- 
rier is bound. This apprpsimation, first developed by 
Ryabchenko and Semenovc^l, allows for the exact calcu- 



lation of the single polaron partition function and makes 
the polaron-pair case much more tractable. 

To consider interactions between two polarons, we 
must make the additional conjecture that there is an in- 
terstitial region between the two polarons within which 
the magnetic ions interact significantly with both carri- 
ers. Such a region must exist in order for the indirect 
ferromagnetic carrier-ion-carrier interactions to be sig- 
nificant. In order to treat the effects of these interstitial 
ions within our model, we assume a cylindrically sym- 
metric interstitial region within which all of the magnetic 
ions interact with both of the carriers with an exchange 
constant, K' . In this interstitial region approximations, 
carrier-ion exchange causes both carrier spins to anti- 
align with the interstitial spins and thereby align with 
each other. Thus, an indirect source of carrier-carrier 
ferromagnetism is introduced into the model. 

In the end, the above approximations yield the model 
Hamiltonian 



if™ = if [(si- Si) + (82-82)] 

-|-ii:'(si-f S2)-S3 + JSi-S2 



(2.4a) 



Si — ^ Si S2 = ^ Sj S3 = ^ Sk 



Interstitial 



(2.4b) 

where K is the intra-polaron ion-carrier exchange con- 
stant, K' is the interstitial ion-carrier exchange constant, 
J is the direct carrier-carrier exchange constant, Si and 
S2 are the carrier spins. Si and S2 are the net polaron 
spins, and S3 is the collective spin of the interstitial re- 
gion. (At this point, we specify only that the interstitial 
region have cylindrical symmetry and be located between 
the polarons. However, for computational purposes, a 
particular shape must be chosen. This matter is dis- 
cussed further in Sec. IV B.) As is indicated in Fig. ^ 
where the details of this model are presented graphically. 
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the essence of the polaron-polaron model reduces to a 
competition between the direct antiferromagnctic carrier- 
carrier interactions characterized by J, and the indirect 
ferromagnetic carrier-ion-carrier interactions character- 
ized by K' . By showing that there are circumstances in 
which the ferromagnetic interaction dominates, a theo- 
retical justification for DMS ferromagnetism can be ob- 
tained. By applying the RS (Ryabchenko and Semenov) 
single step approximation and the interstitial approxi- 
mation to the polaron-pair model, we have effectively 
separated the polaron-pair system into one mechanism 
for polaron formation and another for polaron-polaron 
interaction. It is this separation that makes possible a 
calculation of the polaron-pair partition function. 



III. PARTITION FUNCTION CALCULATION 

To obtain the partition function for this problem, 
we must consider both the interacting part of the sys- 
tem (the carriers and magnetic ions within the polarons 
and the interstitial region), and the non- interacting part 
(those magnetic ions that arc external to both the po- 
larons and the interstitial region). Thus, the full parti- 
tion function takes the form, Z = ZppZg^t, where Zpp is 
the polaron-pair partition function and ZexL is the parti- 
tion function of the non-interacting external ion spins. 

Since the external ions are inert, Z^^t is just given 
by the degeneracy of the magnetic ion spins. Taking the 
magnetic ions to be spin-5/2 (as is the case for Mn), each 
spin has 2s-|-l = 6 orientations. Therefore, Zext = 6^"=°=*, 
where N^xt , the number of external magnetic ions in the 
system, is equal to the total number of ions in the system 
minus the number of ions in the polarons and interstitial 
region. Note that while these non-interacting external 
spins contribute no energy to the system, they contribute 
nonzero entropy and therefore cannot be neglected. 

Due to the approximations made in developing our 
model, the polaron-pair system separates into an individ- 
ual polaron part and a polaron-polaron interaction part. 
Hence, as will be explicitly shown, the polaron-pair par- 
tition function can be expressed as the product of two 
single polaron partition functions and a polaron-polaron 
interaction partition function. 



A. Single Polaron Partition Function 

The single polaron partition function can be calculated 
exactly for the single step model that we have adopted. 
For a single polaron, the Hamiltonian is 



H = X,Ks ■ S = [(s + S)2 - - S^] 



(3.1) 



where s is the carrier spin, S is the sum of the magnetic 
ion spins within the polaron, and As is a placeholder con- 
stant that has been inserted for notational convenience 



and will eventually be set equal to one. For a given S, 
the total spin can take two values: S + ^ or S — The 
former yields an energy and degeneracy 



XsKS 



g+ = 2{S + 1) 



while the latter yields 

XsK{S + l) 



9- = 25. 



(3.2) 



(3.3) 



Thus, the single polaron partition function is given by 

Zp^i = Tr [e-/^^] = Y,D{S) [g+e-^^^ + g.e-^^-] 
s 

(3.4) 

where -0(5*) is the number of ways in which the ion spins 
can be arranged to give a collective spin S. Defining 
Dz{S) to be the number of ways that the ions can be 
arranged to give a collective z-component of spin equal 
to S and doing a bit of algebraic manipulation, we obtain 



57Vi/2 

^ £>,(S)cosh(75) 

S=-5JVi/2 



(3.5) 



where 7 = Xsl3K/2 and Ni is the number of ion spins 
within the polaron. Using the definition of a delta- 
function, Dz{S) can be written as 



D,{S) = Tr5\ S-y^S^ 



dX 
2^' 



[QF{iX)]' 



(3.6) 



where 



F(.) = - 



(3.7) 



For large iVi, the sum over S can be converted to an 
integral and the partition function can be written as 



Zpol 



1 + (e^ - 1) 



J -00^ 



d_ 

97 



dS 



e^^ + e 



Continuing to the imaginary temperature axis, using the 
definition of a delta-function, and continuing back to the 
real temperature axis, we obtain 



Zpol = 6^^Zi Zi 



1 + (e^ - 1) 



d_ 



F{lf' (3.9) 
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where we have separated out the factor of which will 
be cancelled by part of Zg^^t in the full partition func- 
tion. Note that this expression has the correct infinite 
temperature limit since Zpoi{T — > oo) — > 2(6)^^, which 
is the partition function for a non-interacting system of 
A^i spin- 5/ 2 magnetic ions and one spin- 1/2 carrier. 



B. Polaron-Pair Partition Function 

The exact quantum mechanical calculation of the 
polaron-pair partition function is significantly more com- 
plicated than the single polaron case. However, at low 
temperatures, a semiclassical technique introduced in 
Ref. |l^ can be used to find Zpp throughout the temper- 
ature range of interest. Specifically, we must make the 
assumption that the temperature is low enough such that 
the ion spins within the polarons and interstitial region 
are well enough aligned that the two polaron spins. Si 
and 5*2, and the interstitial region spin, 6*3, are large 
enough to be treated as classical magnetic moments. 
Thus, we make a semiclassical approximation in which 
Si, 5*2, and S3 are treated as classical spins while the 
carrier spins, si and S2, are treated quantum mechani- 
cally. For the case we are interested in, appropriate for 
light to moderately doped II- VI based DMS, we expect 
that K' , J <^ K . Consequently, we will first find the par- 
tition function for the case of non-interacting polarons 
{K' = J — 0) and then include the effects of nonzero K' 
and J as first order perturbations. Separating the Hamil- 
tonian into three parts, we write H = H0+H1+H2 where 

i/o = A-[A,i(si.Si) + A,2(s2-S2)] 
Hi = J(si • S2) 

H2 = K' [A,i(si • S3) + A,2(S2 • S3)] (3.10) 

and we have introduced four new constants, Asi, \s2^ Ad, 
and Aci . While these constants will be set equal to unity 
at the end of our calculation, they act as placeholders 
which will be useful when we optimize our model param- 
eters in Sec. IV. 

In the non-interacting polarons limit {K' = J 0), 
the polaron-pair partition function is simply Zpp = 
QN1+N2 ^^^^ where Zi and Z 2 ar e the single polaron par- 
tition functions given by Eq. 3.9 with 71 = XsiPK/2 and 
72 = Xs2f3K/2 respectively. In the semiclassical limit, 
the wave functions of the non-interacting polaron-pair 
are: a(l)Q!(2), a(l)/3(2), /3(l)a(2), and /3(l)/3(2), where 
a is the total spin 5-1-1/2 state of the single polaron 
(carrier and magnetic ions aligned) and /3 is the 5* — 1/2 
state (carrier and magnetic ions anti-aligned) . At the low 
temperatures for which the semiclassical approximation 
is valid, only the ground state, /3(l)/3(2), is significantly 
occupied. Hence, the K' and J perturbations are taken 
to be the diagonal matrix elements in this ground state. 

Therefore, for interacting polarons with nonzero K' 
and J, we write 



Mj = m)m\Hi\mm)) ^i^'-^' ^^^^ 



4 ^1^2 



= {mm\H2\mm)) 



"2" 



AclS'2Si -|- Ac2>S'iS2 

S'lls'2 



•s. 



where 



V2 



c2 



+ AclAc2Ail2, 



4 

(3.11) 



M3S3 
(3.12) 



(3.13) 



/ii2 is the cosine of the angle between Si and iS'2, and 
/i3 is the cosine of the angle between S3 and the z-axis. 
Making use of these matrix elements, the polaron-pair 
partition function can be written as 



Tr 



-(3(Ho + Mj+Mk,) 



X D{Si)DiS2)D{S3)e 



d^Sid^S2d^S3 



7lSl+72S2-^+n/i3S3 



(3.14) 



Performing the indicated integration and doing a bit of 
algebra this becomes 



Zpp = 6^^+^' Z1Z2 



v/2 



e [6F{fl)f''xdx 



(3.15) 

where = 1 + /ii2, N3 is the number of magnetic ion 
spins within the interstitial region, and we have identi- 
fied Zi and Z2 as the single polaron partition functions. 
Notice that in this approximation, the partition func- 
tion does separate into a polaron formation factor and 
a polaron-polaron interaction factor. Multiplying this 
result by the partition function for external ion spins, 
Zext = Q^tot-^i-^^-Na ^ ^^^^ dropping the constant fac- 
tor of 6^'°* , we obtain the full partition function 



Z — Z1Z2Z3 



Z3 = 



V2 



(3.16) 



where Z3 is the polaron-polaron interaction part. 



IV. VARIATIONAL OPTIMIZATION OF MODEL 
PARAMETERS 

Given the partition function calculated in the previ- 
ous section, we proceed to optimize the parameters of 
our model via a variational approach. At T = 0, opti- 
mal values of the model parameters could be obtained by 
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minimizing the expectation value of the model Hamilto- 
nian. At T 7^ 0, we adopt anpapalogous variational ap- 
proach described by Feynmanc^ for which the quantity 
to be minimized is the JF-function 



T — F„i + {H — Hjn) 



(4.1) 



where Hm is the model Hamiltonian and H is the true 
hamiltonian. The aveeage (. . .) is taken over the states of 
Hm- P can be shownnJ to be an upper bound on the true 
free energy F of the Hamiltonian H at the temperature 
T in question. By minimizing T with respect to the 
model parameters, optimal values can be determined as 
a function of temperature. 



As will be shown explicitly in Sec. IV B , the total T- 
function separates into the sum of two single polaron 
functions and a polaron-polaron interaction function. 
Therefore, we shall optimize the single polaron param- 
eters first and then consider the interaction parameters. 



A. Single Polaron Parameter Optimization 

The single polaron J^-function can be obtained by 
expressing (-ffm), (^), and Fm in terms of the single 
polaron partition function, Z\. Recall that the model 
Hamiltonian has the form 



Fir, 



(4.2) 



where index j runs over all magnetic ion spins, Sj , within 
a sphere of radius R about the carrier spin, s, and As is 
a constant which will soon be set equal to unity. Taking 
the thermal average over the eigenstates of Hm yields 
that 



{Hm) = XsKNi{s-S.)=- 



Xs dlnZ-i 
J dXs 



(4.3) 



where Ni is the number of magnetic ions within the po- 
laron. The true Hamiltonian has the form 



K'Vol 



0.001 



R vs.T 





FIG. 4. Temperature dependence of single polaron param- 
eters. All distances are in units of a_B and all energies are in 
units of a/a% ~ 625 K. 



(4.4) 



where the index n runs over all magnetic ion spins. 
Therefore, noting that (s • S„) is only nonzero for ion 
spins within the polaron, we find that 



(H) 



7s dlnZi 



where 



Is 



(4.5) 



(4.6) 



and the integral is over a sphere of radius R and volume 
Vs- Finally, since the free energy is 



Fm = --InZi 

P 



^1 



1 
7 



lnZi-|-(7s-l) 



SlnZi 

dXs 



(4.7) 



we can combine the expressions above (setting As — 1) 
to obtain the single polaron JF-function 



(4.8) 



where Zi is given by Eq. ( p.9| ). By minimizing with 
respect to the parameters R and K , the optimal values 
of these parameters can be found. 

Performing such a procedure numerically over a range 
of temperature values, the optimal values of the single 
polaron parameters were determined as functions of tem- 
perature. The results of this optimization for a magnetic 
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ion density of 5 ions per sphere of radius las are plotted 
in Fig. 

In the high temperature hmit (T 3> K), the exchange 
interaction between the carrier and the magnetic ions 
within the polaron is insignificant compared to temper- 
ature. Thus, the magnetic ion spins are not ahgned and 
there is no difference between the free energy of a spin 
within the polaron and that of an external spin. As a re- 
sult, the J-"-function is minimized when the model carrier 
wave function best matches the true carrier wave func- 
tion. This matching of a step of width R and height 
K to a hydrogenic wave function yields the optimal, 
temperature-independent values of R and K . Thus, as 
is shown in Fig. |[ i?, K, and KVs (the total exchange 
energy) are temperature-independent in the high temper- 
ature regime. 

For low temperatures {T <C K), the carrier- ion ex- 
change interaction is significant compared to tempera- 
ture. Thus, the magnetic ion spins located near the car- 
rier anti-align with the carrier spin. In this situation, the 
inclusion of an additional ion within the polaron entails 
a gain in exchange energy. However, since the number of 
external spins decreases by one, there is also a decrease 
in the entropy of free spins. Therefore, the optimal R is 
determined by the balance of exchange energy and the 
entropy of free spins. As T decreases, the exchange en- 
ergy gained by increasing the size of the polaron becomes 
more valuable. Thus, as is shown in Fig. ^, R increases 
as log(l/T) as T drops. As R increases, K must decrease 
in order to maintain the match between the model step 
wave function and the true wave function. Thus, K de- 
creases with decreasing T. Finally, despite the decrease 
in the total exchange energy (which is proportional 
to KVs) increases as T drops and the spins align. 

It is interesting to note that in this problem, the vari- 
ational principle leads one to match the Hamiltonian at 
high T, while entropy-energy balance determines the pa- 
rameters at low T. This is just the converse of what one 
expects in most problems. 



B. Polaron-Pair Parameter Optimization 

Just as for the single polaron, we can obtain the 
polaron-pair J^-function by expressing (i?™), {H), and 
Fm in terms of the full partition function, Z . For the 
polaron-pair, our model Hamiltonian is 

j 

+ XclK' ^ Si • Sfc + \c2K' S2 • Sfc + JSi • S2 

k k 

(4.9) 

where indices i and j run over the magnetic ions in po- 
laron 1 and 2 respectively, k runs over ions in the cylin- 
drically symmetric interstitial region, and the A's are con- 



stants which will soon be set equal to one. Taking the 
thermal average then yields 

(i?„,) ^ XsiKNiisi ■ S,) + Xs2KN2{s2 ■ Sj) 
+ X^iK'Nsisi ■ Sk) + Xc2K'N3{s2 ■ S^) 
+ J(si • S2) (4.10) 

where iVi = N2 are the number of ions in polaron 1 and 
2, and A'3 is the number of ions in the interstitial region. 
Noting that 



(Si • S,> - (S2 • S,) 
(Si • Sfe) = (S2 • Sk) 



1 IdlnZ 



KN, p dXsi 
1 ld\nZ 



and setting the A's equal to one, this becomes 



(4.11) 



d\nZ d\nZ 



dXsi dXci 
Since the true Hamiltonian has the form 



H = a}_^si- S„|^'i„|^ + a2^S2 • S„|*2n| 



J(si-S2). (4.12) 



JSi • S2 



(4.13) 



and we know that (s • S„) is only nonzero for ion spins 
within the polarons or interstitial region 

{H) = a(si • S,) I^HP + • Sfe) Y 1*1-^1" 

i k 

+ a{s2 ■ S,) Y l*2,f + a{s2 ■ Sk) ^ |*2feP 



J(si • S2). 



(4.14) 



Again making use of Eq. (4.11), this becomes 



(H) 



Is 



d\nZ 
dXsi 



+ 7c 



d\nZ 



dXci 



J(si-S2) (4.15) 



where 7s is defined in Eq. (4.6) and we have now defined 

(4.16) 



7c 



— f d\\^i(r) 



where the integral is over the cylindrically symmetric in- 
terstitial region of volume Vc- Subtracting Eq. ( 4.12| ) 
from Eq. ( 4.15 ) and adding the free energy, Fm = 
— ^ InZ, we obtain the polaron-pair J^-fimction 



.d\nZ , , In Z 
lnZ + 2(7,-1)-— + 2(7,-1) 



dXs\ 



dXci 
(4.17) 



where Z is the partition function given by Eq. ( ^.16 ). 
Since Z is of the form Z = Z1Z2Z3, it is clear that J-" is 
of the form = !Fi + !F2 + ^3 and therefore separates 



7 




FIG. 5. Schematic of ellipsoidal lens shaped interstitial re- 
gion. 



into a polaron formation term, J-2 — 2 J^i , which we 
developed in the previous section, and a polaron-polaron 
interaction term 



lnZ3 + 2(7e-l) 



91nZ3 

9A, 



cl 



(4.18) 



Since these terms share no variational parameters, the 
polaron-pair model can be optimized by using Ti to op- 
timize K and R (as we did in the previous section) and 
using to optimize K' and the parameters describing 
the geometry of the interstitial region. 

Before we can proceed to minimize J-3, we must define 
a specific geometry for the interstitial region between the 
two polarons. The bi-spherical geometry of the problem 
suggests that a natural choice would be the spherical lens 
formed by the intersection of two spheres centered on the 
two polarons. However, such a shape can be completely 
specified by a single parameter, the lens width h. In order 
to provide an additional degree of freedom within the 
model, we will use the slightly more general ellipsoidal 
lens formed by the intersection of two ellipsoids centered 
on the polarons. In this manner, the interstitial region 
can be specified by two parameters, the lens width, /i, and 
the lens radius, p. The details of this shape are depicted 
in Fig. I 

The task of optimizing the model parameters, K' , h, 
and p, is complicated in two ways. First of all, unlike the 
single polaron case where both the partition function, Zi, 
and the geometrical factor, 7^, could be obtained analyt- 
ically, both the interaction partition function, Z3, and 
the interstitial geometrical factor, 7c, must be calculated 



numerically. This complicates the numerics but poses no 
fundamental problem. 

The second complication requires a bit more atten- 
tion. Naively one would expect that by blindly vary- 
ing these three parameters until J^^ is minimized, the 
optimal values of these parameters could be obtained. 
However, upon closer examination, it becomes clear that 
this is not the case. In performing this optimization of 
the interstitial region parameters, it is our objective to 
determine the parameters that best match the true in- 
teraction between both carriers and the magnetic ions 
in the interstitial region. However, minimizing the JF- 
function merely yields the configuration that, overall, is 
most energetically favorable. As a result, if all three pa- 
rameters are varied, the minimum J^^ will be achieved 
when the interstitial region has totally engulfed both of 
the polarons. In this configuration, the interstitial re- 
gion contains ions that are very close to the carriers and 
therefore experience large exchange interactions. How- 
ever, these ions that interact very strongly with one of the 
polarons barely interact at all with the other. Thus, the 
result of an unconstrained minimization does not yield 
an interstitial region with which both of the carriers in- 
teract strongly. Therefore, to obtain sensible results, the 
minimization of the interaction ^-function must be con- 
strained. Two methods of constraining the minimization 
have been developed: the fixed width method and the 
optimally spherical method. 

The interstitial region in which both carrier wave 
functions are significant must be concentrated about 
the point halfway between the two polarons. Thus, 
one method of constraining the minimization of the J-- 
function is to fix the lens width, h, to a set value. Using 
this fixed width method, the only parameters that are 
allowed to vary are the interstitial interaction strength, 
K', and the lens radius, p. For a given temperature, the 
jF-function will be minimized for some optimal values of 
these two parameters. Thus, by performing a numeri- 
cal minimization for a range of temperature values, the 
temperature dependence of K' and p can be determined. 

In the fixed width model, the fixed value of h is chosen 
such that both carrier wave functions will be "significant" 
within the interstitial region. In a sense, the choice of a 
particular h defines the threshold of exchange interac- 
tion strength with both carriers that is required for an 
ion to be included in the interstitial region. However, 
as temperature changes, this threshold should change as 
well since all energies in the system are measured with 
respect to temperature. Thus, as temperature drops, the 
threshold should also drop and the width of the inter- 
stitial region should increase (to include those additional 
ions that now meet the lowered threshold) . This is a fea- 
ture that is not incorporated into the fixed width method, 
since h remains constant as temperature varies. Unfor- 
tunately, this problem cannot be solved by making h a 
variational parameter since, if this is the case, the min- 
imum JF-function will only be achieved when the inter- 
stitial region engulfs the polarons. One way to incorpo- 
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rate the growth of h with decreasing T into the process 
would be to choose a new fixed h for every T and set the 
temperature-dependent threshold for inclusion in the in- 
terstitial region by hand. However, this approach would 
require a prior knowledge of the temperature dependence 
of h, which we do not have. 

A more natural way to allow the threshold to vary 
with T is to let the threshold be set by the geometry 
of the system. The polaron-pair problem consists of two 
sources of radially symmetric wave functions separated 
by an inter-impurity distance. Thus, the system has an 
innate bi-spherical geometry. If one were to define a re- 
gion in space in which the wave functions of both car- 
riers would surpass a given threshold, the geometry of 
the system and the spherically symmetric nature of the 
wave functions would dictate that this region be the in- 
tersection of two spheres centered at the carrier sites - a 
spherical lens. Therefore, the natural shape for the in- 
terstitial region is a spherical lens. This fact provides the 
condition for setting the width of the interstitial region 
in what we call the optimally spherical method. For a 
given temperature, the fixed width, h, is set to be that 
for which the minimization of the J^-function automat- 
ically yields a value of the lens radius, p, for which the 
ellipsoidal lens becomes spherical. Thus, the result is a 
spherical lens interstitial region in which the lens radius is 
the energetically optimal radius for the given width. As 
we will show, this technique yields the proper increase in 
h as temperature drops. 

Using both the fixed width approach (with h set to 
2aB) and the optimally spherical approach, the polaron- 
polaron interaction parameters, K' , h, and p, were opti- 
mized as functions of temperature. The results of this 
optimization for a magnetic ion density of 5 ions per 
sphere of radius las and a carrier density such that 
Ri2 — 6a B are plotted in Fig. ^. Although the parame- 
ter values are plotted over a wide range of temperatures, 
they are only meaningful for temperatures at which the 
polaron-pair model is valid. Clearly, once the polarons 
have increased to a size such that they touch the inter- 
stitial region, the polaron-pair model is no longer valid. 
This impact between the polarons and the interstitial re- 
gion occurs at T = Timpact where 2R + h = i?i2. In 
Fig. ^, the temperature of impact has been denoted by 
a vertical dashed line for both the fixed width and the 
optimally spherical cases. Note that for either method, 
^Impact is found to be approximately equal to 0.005 a/a^ 
or around 3 K for typical ion and carrier densities. 
Since K' (Ti^npact) ^ ^Impact, this cutoff temperature lies 
within the high temperature regime (T ^ K') where the 
interstitial ion-carrier interactions are insignificant com- 
pared to temperature. Hence, for temperatures where 
our model is well defined, all of the polaron-polaron in- 
teraction parameters are constant. For very low tem- 
peratures (T <C K'), our results are not quantitative, 
but since polarons are nearly aligned (either parallel or 
antiparallel) by this point, this regime lies beyond the 
temperature range of interest. The use of a constant pa- 



■ Fixed Width 

o Optimally Spherical 



FIXED WIDTH 
IMPACT TEMP 



SPHERICAL 
IMPACT TEMP 




FIG. 6. Temperature dependence of polaron-polaron inter- 
action parameters obtained using both the fixed width and 
optimally spherical methods of constraining the minimization. 
All distances are in units of as and all energies are in units 
of a/o| ~ 625 K. 
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rameter model to treat the polaron-polaron interaction 
in Ref. n9^ is therefore justified. 



V. POLARON-POLARON INTERACTION 

Given the polaron-polaron interaction partition func- 
tion, Z3, as expressed in Eq. ( 3.16 ), it follows that the 
thermal average of the cosine of the angle between the 
two polaron spins can be obtained via 



/c 



4kgT p 



(COS 6*12) = 



K'x 
V2kBT 



AT, 



X dx 



rV2 - 

Jo e 



J(x 



'F 



( K'x \ 
\V2kBT J 



N3 



X dx 



(5.1) 



where F{z) is defined in Eq. (3/7), — 1 = cos6'i2, and 
7V3 = Ns^h, p) is the number of magnetic ions in the inter- 
stitial region. As motivated above, the polaron-polaron 
interaction parameters, K' and A^3, are approximately 
temperature-independent in the temperature range of in- 
terest. Thus, for various constant values of J, K\ and 
A^3, the above expression can be evaluated numerically 
as a function of temperature to reveal the temperature 
dependence of the inter-polaron angle. The results have 
been plotted in Fig. |^ for A^3 = 20 and several values 
of the ratio K' / J . We see that for reasonable parameter 
values, (cos 612) is positive and neighboring polarons tend 
to align. For large enough K' / J, indirect ferromagnetic 
carrier-ion-carrier interactions can dominate the direct 
antiferromagnetic carrier-carrier interactions to yield a 
net ferromagnetic polaron-polaron interaction. 

Polaron Pair - Constant Parameters - K7J={0,0.05,0.10,0.15} 




VI. CONCLUSIONS 

In this analysis of the magnetic behavior of diluted 
magnetic semiconductors, we have proposed a simplified 
model to describe both the formation of bound magnetic 
polarons and the interactions between them. Approxi- 
mating carrier wave functions via sharp cutoffs that de- 
fine polarons and interstitial regions between them, we 
have obtained a tractable model and calculated the re- 
sulting partition function. Utilizing a finite temperature 
variational approach, the model parameters have been 
optimized as functions of temperature. At very high 
temperatures (T ^ K), the spins within the system are 
not aligned and the model parameters obtain the con- 
stant values for which the model wave functions best 
match the true wave functions. At lower temperatures 
{K' <C T <C K), where the carrier- ion exchange inter- 
action becomes more significant, both polaron size and 
total intra-polaron exchange energy increase logarithmi- 
cally with decreasing T. However, throughout this inter- 
mediate temperature range, the parameters controlling 
polaron-polaron interactions, K' , h, and p, remain con- 
stant as T varies. At very low temperatures (T ^ K')^ 
even these interstitial region parameters would become 
temperature dependent. However, for reasonable car- 
rier and ion densities, the polarons grow large enough to 
touch the interstitial region, such that our model ceases 
to be valid, well before this regime is realized. Hence, for 
temperatures of interest, a model in which the polaron- 
polaron interaction parameters are constant is justified 
variationally. Therefore, for reasonable values of K' / J 
we obtain, as in Ref. 09, a net ferromagnetic polaron- 
polaron interaction, in agreement with the experimental 
evidence for ferromagnetism in insulating doped DMS. 

The further growth of the polaron pair bubble, for 
T <C K' ^ could be captured via a variational scheme em- 
ploying an appropriately generalized object. Starting at 
the point of polaron overlap, such an object, involving 
two or three variational parameters, should evolve from 
a "peanut shape" , with cylindrical symmetry, to a sphere 
with the mid-point of the two dopant sites as its center. 
This could be a promising direction for future research. 

Since this work was completed, a tendency towaxdJer- 
romagnetism has been observed by other groups.L3c3 In 
addition, a separate mechanism for ferromagnetic align- 
ment, resulting from the strong local exchange fields ex- 
perienced at the two impurity sites due to nmrby Mn, 
has been considered by Angelescu and Bhatt.EJ 



FIG. 7. Thermal average of the cosine of the angle be- 
tween polaron spins, (cos 612), plotted versus temperature for 
A^3 = 20 and several values of the ratio K' / J. 
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